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It is understood that congestion in traffic can be interpreted in terms of the in¬ 
stability of the equation of dynamic motion. The evolution of a traffic system from 
an unstable or metastable state to a globally stable state bears a strong resem¬ 
blance to the phase transition in thermodynamics. In this work, we explore the 
underlying physics of the traffic system, by examining closely the physical properties 
and mathematical constraints of the phase transitions therein. By using a meso¬ 
scopic approach, one entitles the catastrophe model the same physical content as 
in the Landau’s theory, and uncovers its close connections to the instability of the 
equation of motion and to the transition between different traffic states. In addi¬ 
tion to the one-dimensional configuration space, we generalize our discussions to the 
higher-dimensional case, where the observed temporal oscillation in traffic flow data 
is attributed to the curl of a vector field. We exhibit that our model can reproduce 
the main features of the observed fundamental diagram including the inverse-A shape 
and the wide scattering of congested traffic data. When properly parameterized, the 
main feature of the data can be reproduced reasonably well either in terms of the 
oscillating congested traffic or in terms of the synchronized flow. 
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I. INTRODUCTION 

Traffic flow modeling has long attracted the attention of physicists (for reviews, see for 
examples nHEi). In particular, it is understood by many authors that the traffic congestion 
is closely connected to the instability of the equation of motion (EoM) of the corresponding 
traffic theory [ZHn]. For instance, by using a microscopic car-following model, Bando et 
al. [H] investigated the imaginary part of the frequency of small oscillations around a given 
steady state. It turns out that the signature of the imaginary frequency indicates the 
properties of linear stability of the steady state solution: negative imaginary frequency 
implies that small perturbation will blow up exponentially to render the system unstable. 
For hnite size perturbations, one usually has to resort to numerical techniques. By numerical 
calculations, the authors showed that hnite size perturbation will cause the vehicles to 
cluster into congestion, and further study showed m that in this specihc model the traffic 
congestion is related to the continuously deformable kink soliton solution. Another study 
was carried out by Kerner and Konhauser [7] from a hydrodynamical point of view. Cluster 
formation was encountered numerically where small deviations about homogeneous steady 
traffic how destabilize the system. Although the above approaches are diherent in their 
respective levels of aggregation, they are very similar in nature since they all follow the 
general arguments used to investigate the stability in a classical dynamical system na. 

When the system is locally unstable for inhnitesimal oscillations, it is not known be¬ 
forehand whether hnite deviation from a steady state may lead to global instability and 
consequently to any unphysical situation. For a realistic approach, the system must always 
be bounded and therefore, divergence shall not occur to any physical observable. For in¬ 
stance, quantities such as the total number of vehicles, vehicle velocity and acceleration must 
always remain hnite; in other words, the system can be linearly unstable, but the amplitude 
shall never increase unboundedly. In fact, such physical constraints can be seen as necessary 
conditions for a realistic model. If one follows this train of thought, it is noted that for a 
thermodynamical system the phase transition is usually caused by perturbations at locally 
unstable or metastable states. Under such situation, the system evolves and eventually 
relocates itself to another locally or globally stable state; at the same time, the thermo¬ 
dynamical system in question must always be bounded. To determine if one state is more 
favorable than another, usually the thermodynamical potential is introduced, so that the 
locally stable and metastable states are related to the local minima while the globally stable 
state corresponds to the global minimum of the potential. The well-known Landau’s theory 
of phase transition [13] is one of such phenomenological approaches dealing with transitions 
between different phases via the concept of free energy. 

Incidentally, in the traffic flow theory the catastrophe model HI bears an analogy with 
the Landau’s theory. The theory was hrst developed by Thom [15] and Zeeman [16] and was 
used mostly to model catastrophic events. To describe the observed discontinuous transition 
between different traffic states, it was later applied to the held of traffic how. To this end, 
a potential function is introduced, where the minima of the function correspond to the 
observed traffic states. In the case of the cusp catastrophe model HU, the extrema of the 
function are single-valued on one side but three folded on the other. The folded surface 
in the parameter space as well as the consequent occurrence of sudden jump of observables 
featured in this model are indeed very similar to those in the theory of phase transition. 
Despite its success, the catastrophe model does not have a microscopic origin; it attempts 
neither to deal with the temporal evolution of the traffic system, nor to derive the form of the 
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potential function from first principle. Naturally, there is no concept such as temperature or 
thermodynamical ensemble in the catastrophe approach, so the physical content of the model 
seems to be very different from that of Landau’s theory of phase transition. Nevertheless, 
the main idea proposed by the catastrophe model is intriguing owing to its unique viewpoint 
to the problem. 

Over the last decade, the investigation of traffic ffow was advanced by more realistic ap¬ 
proaches featuring parametrization based on observed traffic data, vehicle models as well as 
government regulations. Many ingenious works [3 El [IHH22] were carried out aiming at a 
deeper insight into the traffic system; in particularly, the study of the different phases of the 
traffic system and their coexistence and transitions. As an example, in [22], a microscopic 
intelligent driver model (IDM) was proposed and used to explore various phases of traffic 
ffow. By numerical simulations, the model was shown to be very successful in reproduc¬ 
ing the observed data once the boundary conditions of the simulated system were properly 
given. For a homogeneous system, both the free traffic (FT) state and the homogeneous 
congestion traffic (HCT) state were reproduced by the model and are shown to be linearly 
stable against small perturbations. What makes the model more signihcant is when it is 
applied to study traffic systems with inhomogeneity. The authors show that the model is 
capable of reproducing other non-st at ionary traffic states, in particular, the oscillatory states 
observed in the data [23fl26] which are characterized by prolonged temporal oscillation of 
physical quantities, such as oscillating congested traffic (OST) and triggered stop-and-go 
waves (TSG). In the literature, such phenomena were also studied and reproduced by other 
studies [HEllIHHSP. In fact, the periodic orbit in dynamical system is a topic of increasing 
interest in various areas of science, known as limit cycle [T2|. In the case of traffic ffow, 
they are understood by most authors as distinguished traffic states m It is worth noting 
that in the context of traffic ffow, the physical content of the above oscillatory phases is 
usually attributed to the inhomogeneity of the traffic system [H EH] and are treated by in¬ 
troducing either some extra external source (such as on-ramp flow) into the system, or local 
discontinuity into the model parameters. As far as the authors are aware, in the literature 
most studies were for specihc traffic scenarios with appropriately chosen parameters as well 
as initial/boundary conditions. Besides, as most results were drawn based on numerical 
simulations with properly chosen initial and/or boundary conditions, the properties and 
interpretation of the resulting traffic state sometimes vary in different models. In addition, 
although the temporal oscillatory traffic state is generally recognized as a traffic phase, its 
connection to the concept of phase transition in thermodynamics is somewhat obscure. This 
motivated us to elaborate the starting point of our approach of the traffic system from the 
general concepts of stability and of phase transition, rather than from any particular charac¬ 
teristic of the traffic system. The non-stationary nature connected to the time dependence 
of measured physical quantities seemingly rules out the possibility to employ a theoretical 
framework which only deals with quasi-static process and/or steady state, such as the the¬ 
ory of phase transition in thermodynamics. Nonetheless, it remains unclear and therefore 
captivating whether such oscillatory states can also be understood by an approach following 
the very spirit of the potential function. 

The present work aims at a dynamical traffic theory, which reproduces the main char¬ 
acteristic of observed traffic flow data but expressed in the language of the catastrophe 
model. In particular, we study stationary as well as non-stationary traffic phases and the 
transitions among these states by investigating the stability of the corresponding EoM of 
the traffic system. In our model, the EoM is based on the kinetic theory of gas, namely. 
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a mesoscopic Boltzmann-type theory. It is shown that the potential function can be de- 
hned in terms of the transition coefficients of the Boltzmann-type equation. For the case 
of steady state traffic, the process to minimize the potential function is equivalent to pick 
out the stable steady states of the EoM. In this context, our approach “connects the dots” 
between the catastrophe model and Landau’s theory of phase transition from the viewpoint 
of a mesoscopic dynamical theory. Moreover, our theory explicitly depends on time, which 
makes it feasible to investigate the temporal evolution of system and the possibility of any 
(stable) oscillatory solution. For the most part, our study strives to establish a general 
framework of a dynamical theory rather than to focus on specific parameterization. The 
paper is organized as follows. In section II, we give a brief review on the traffic flow theories, 
and then we present the mesoscopic traffic flow model. In section III, we investigate steady 
state traffic, the properties of instability of the EoM, as well as realistic physical constraints. 
We show that the present model possesses most features of the catastrophe model once the 
potential function is properly defined in terms of the transition coefficients. Furthermore, 
the potential function of the model can be seen on the same footing with the thermody¬ 
namical potential and the relation to the Landau’s theory of phase transition is established. 
The model in one-dimensional configuration space is investigated. In section IV, in order 
to study non-stationary states, we generalize the model to the higher-dimensional case. It 
is shown that such complication provides the possibility to model the observed temporal 
oscillation of traffic flow states. We investigate two possible scenarios: in the first scenario, 
the observed oscillation in traffic flow data is understood to be related to the limit circle 
solution in configuration space; in the second scenario, the potential function is parameter¬ 
ized to possess a metastable region in accordance to the definition of the synchronized flow 
introduced in the three phase traffic theory [21128]. In both cases, the temporal oscillation 
is attributed to the curl of a vector potential. The fundamental diagram in both cases are 
also studied. The last section is devoted to discussions and concluding remarks. 


II. A MESOSCOPIC TRAFFIC FLOW MODEL 

One important empirical measurement for a long homogeneous freeway system is the 
so called fundamental diagram of traffic flow. It is plotted in terms of vehicle flow as a 
function of vehicle density (or concentration). The flow-concentration relation possesses the 
following main characteristics: (1) It is divided into two different regions of lower and higher 
vehicle density respectively, which correspond to free and congested flow; (2) The maximum 
of the free flow occurs at the discontinuous junction above that of the congested flow region 
which forms an inverse-A shape and (3) Congested flow presents a broader scattering of the 
data points on the flow-concentration plane. Recently, we proposed a mesoscopic model 
for the traffic flow [29| by the method of stochastic differential equation (SDE) [30], where 
the dynamics of the system is governed by the temporal evolution of the distributions of 
vehicles among different speed states via transitions in configuration space. In addition to 
the conventional transition terms, stochastic transition is introduced in the model to describe 
the stochastic nature of traffic flow. It was shown that analytic solutions can be obtained 
not only for the expected value of speed and traffic flow, but also for their variances. The 
latter is understood to account for the resulting scatter of the data point of the fundamental 
diagram. In particular, we show that a simplified version of our model by assuming only 
two speed states with constant transition coefficients is able to reproduce the empirical 
data on a Brazilian highway [31], and two of the above characteristics, (1) and (3), were 
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well-reproduced. 

In the following, we briefly review our traffic model, more details of the model can be 
found in [29]. Let us consider a section of highway where the spatial variation of the ve¬ 
hicle distribution can be ignored. For simplicity, only discrete values for speed, namely 
^ 1 ,^ 2 , • • • ,vd are considered, and the number density of vehicles conducting at speed state 
Vi is denoted by n^. In time, a vehicle with speed Vi may transit to another state Vj according 
to the following set of SDE 


drii 

dt 


D D 

i=i t=i 


( 1 ) 


where the speed transition on the r.h.s. of the equation is a summation of two contributions: 
the deterministic and stochastic transitions measured by the transition rates Cij and Sij 
respectively. For the stochastic transition, randomness is introduced through the white noise 
Wj, which is a random signal characterized by a featureless power spectral density. When 
j ^ i, the coefficients Cij and measure the rate for a vehicle with speed Vj to transit to 
another state with speed Vi. It has been shown [29] that these transition coefficients are not 
completely independent, and the number of degrees of freedom of the system is D — 1. 

It is noted when one is only interested in the temporal evolution of the expected value of 
physical quantities, as in many cases of physical interest, the stochastic terms can be ignored 
[30] . By removing the stochastic transition terms, Eq.Q reduces to the corresponding 
deterministic EoM 


drii 

dt 


D 


( 2 ) 


As discussed above, different traffic ffow phases are originated from the properties of stability 
of Eq.Q. For the present model, in principle, there are at least two sources of instability. 
The first one is due to the stability of the deterministic EoM Eq.([^, which is completely 
determined by the deterministic transition rate c^. The second source is due to the stochastic 
fluctuation, measured by the white noise Wj and its strength in terms of the transition rate 
Sij. In this work, we will concentrate on the first source. The effect of the stochastic 
transitions was studied in [29] . 

To discuss the stability of the deterministic EoM Eq.([^, one can introduce a small 
perturbation around a steady solution. In general, the transition coefficients are not expected 
to be constants. Nevertheless, for a linear system, they usually can be expanded around 
those of the steady solution. Intuitively, for any physical solution, Ui must be bounded 
from above and below. As a result, a stable steady solution implies that the transition 
matrix, the {D — 1) x {D — 1) matrix defined by the transition coefficients, must be locally 
negative dehnite. This is because, any positive eigenvalue would imply that the vehicle 
number of some state increases exponentially in time near the vicinity of the steady solution, 
which amplihes any small perturbation. On the other hand, for a realistic traffic model 
which provides a natural description of traffic congestion, its EoM must contain region with 
unstable solutions. This is the main topic to be explored in the subsequent sections. 
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III. POTENTIAL FUNCTION IN A ONE-DIMENSIONAL BOLTZMANN-TYPE 

TRAFFIC FLOW THEORY 


In this section, let us first concentrate on the simplest case by employing a simplified 
version of the model which contains only two speed states: Vi and V 2 with corresponding 
occupation densities ni and n 2 - Using this simple model, our goal is to elaborate, to the 
greatest extent possible, a traffic flow theory which consists of different traffic phases and 
the transition between them based on the properties of the instability of the corresponding 
EoM. The equation of motion of this simplified model reads 


drii 

dt 

dn2 

dt 


-Piini + pi2n2 
-P22n2 + P2ini 


(3) 


where one assumes vi < V 2 without losing generality. When adopting properly chosen 
parameters [29], the above simple model has be shown to reproduce the observed funda¬ 
mental diagram, as well as its observed variance once the noise terms are present. By 
considering the conservation of total vehicle numbers, one has ni + n 2 = N, and therefore 
Pii = -cii = P 21 = C 21 and P 22 = -C 22 = P 12 = C 12 {cij of Eq.Q). The EoM is in fact of one 
dimension: 


dni 

dt 


-{pii +Pi2)ni +P 12 N 


(4) 


For steady solution n 


( 0 ) 

1 ) 


one has 


n 


( 0 ) 

1 


P12 

Pii + P12 


N 


(5) 


For a small perturbation around the steady solution, one can expand it as ni = + 

When the transition coefficients are simply constants, the resulting EoM of small 
perturbation reads 


=-{Pu+Pi2)n^i^ ( 6 ) 

from which one immediately sees that the small perturbations die out exponentially if pn + 
P 12 are positive. 

However, a model with constant transition coefficients are less physically interesting since 
they represent a system which is always stable or unstable. In a realistic scenario, a traffic 
model shall possess different regions featuring distinct stability properties. Fortunately, 
the above argument can be easily generalized to include those elements into the model by 
considering the the case where transition coefficients are functions of ni. In particular, it 
will turn out to be useful to define a potential function, U = U{ni), in terms of the transition 
coefficients by its derivative U' as follows. 


U\ni) = {pii+pi 2 )ni -P12N 


(7) 


In other words, the potential function is determined, up to an irrelevant constant, as follows: 
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FIG. 1. (Color Online) Schematic potential functions for the one-dimensional case: (a) Stable free 
traffic (FT) state at vehicle density N = 0.1, where the FT state is the only stable steady state; (b) 
Stable FT state and metastable homogeneous congestion traffic (HOT) state at = 0.35, where 
both the FT and the HOT states are locally stable solutions, but the FT state is globally stable; 
(c) Metastable FT state and stable HOT state at = 0.4, similar to the case (b) except that 
the HOT state is globally stable; (d) Stable HOT state at A^ = 0.6, similar to the case (a), where 
the HOT state is the only stable steady state. The analytic expression of the potential function is 
given in the Appendix. 


/ ni 

dn[ {{pu{n[) +Pi 2 {n[))n[ - pi 2 {n[)N) (8) 

The potential function defined above is motivated by the fact that for a steady state, 
one has ^ = 0, which implies U' = 0. We are about to show that a locally stable state 
corresponds to a local minimum of the potential function, therefore the potential function 
defined in Eq.(scalarpot) has the exactly same physical content as that of the catastrophe 
model. First, it is easy to verify that the equation for small perturbation, Eq.([^, can be 
generalized and rewritten in terms of the potential function: 


dt 




-U”n\ 


( 9 ) 


where one has U” = (pn + P 12 ) — derivatives are carried out 

with respect to rii and higher order terms are ignored. From Eq.(|^, one sees that small 




















perturbation dies out when U is convex around the steady state. It implies that a stable 
steady state corresponds to a point with U' = 0 and U" > 0, namely, a local minimum of the 
potential function. Now we are in a position to study the implications of the above potential 
function on traffic ffow theory. Since one expects that the model have different stable states 
for different traffic scenarios, it implies: (1) when the overall vehicle density is small, FT 
state must correspond to a stable state; (2) when the overall vehicle density is large, the 
system is likely to be “trapped” in a HCT state. Different traffic scenarios discussed above 
imply that the potential function can not be hxed but it should depend on some quantities 
that characterize the system, such as the overall vehicle density. In addition, at very low 
vehicle density, FT must be a global as well as unique stable state, as one never spots a traffic 
jam at mid-night. And the same arguments are also valid for the stability of the HCT state 
at very high vehicle density. When a steady state is globally stable, it means that there is a 
unique local minimum in the corresponding potential function U. The potential functions in 
these cases are illustrated in Figj^a) and (d) with different total vehicle density N, where 
N is normalized to a range between 0 to 1. Now we consider an intermediate vehicle density, 
where a steady FT state characterized by bigger traffic ffow and a HCT state characterized 
by smaller traffic ffow coexist. In this scenario, both states are locally stable, but only 
one of them is globally stable. In other words, the system locates in a local minimum of 
the potential function instead of the global minimum, namely, a metastable state. It is 
stable against small perturbation until the magnitude of the perturbation exceeds a certain 
critical value, which is determined by the depth of the local minimum. This scenario is in 
fact well-known and discussed by many authors in terms of different traffic models [201122] • 
In our approach, it is naturally implemented in terms of the potential function, which is 
shown illustratively in Figl.(b) and (c). We note that the coexistence of metastable and 
stable traffic states brings up an intuitive explanation of the discontinuity presented in the 
fundamental diagram as shown below. It is also not difficult to see that the above discussions 
fall in the very same line of arguments as the theory of transition of Landau, the potential 
function can be viewed to play the same role as the thermodynamical free energy. 

Moreover, one has to take into account the physical constraints as boundary conditions. 
Since the potential function U is dehned to be a function of ni, obviously the system shall 
only evolve inside the region where 0 < rii < A^. The latter implies that ^ > 0 when 
rii —)■ 0 and ^ < 0 when ni —>■ N, which consequently apply restrictions on the shape 
of the potential function: U' < 0 when Ui —)■ 0 and U' > 0 when ni —)■ N. In the 
one-dimensional case, this means that the potential function must be conhned between its 
boundaries at ni = 0 and ni = N. In the case when one only considers two traffic flow 
phases, the shape of the potential function is unambiguously determined and it shall look 
like those depicted in Figj^ 

As an example, in Figj^a) and (b), we calculate numerically the temporal evolution of 
the system with different initial conditions. We adopt exactly the potential functions shown 
in Fig0a) and (b). One observes that in the case of Figj^a), the system always evolves into 
the globally stable steady state independent of its initial conditions; in the case of Figj^b), 
the system displays the above mentioned behavior for a metastable state. The specihc form 
of the potential function used in Figjljj^can be found in the Appendix. By identifying the 
form of the potential function, it is made possible to determine the mimimal strength of the 
pertubation in order to knock the traffic flow from a metastable state into a global state, in 
other words, the depth of a local minimum of the potential function provides ns quantitative 
information on the criterion of the phase transition. 
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FIG. 2. (Color Online) Temporal evolution of the system for different potential functions in the 
one-dimensional case: (a) Temporal evolution with different initial conditions for the potential 
function shown in Fig{^a). It is noted that the system always evolves into the globally stable FT 
state, independent of the initial condition; (b) Temporal evolution with different initial conditions 
for the potential function shown in Figj^c). It is noted that the system may evolves into either 
one of the two locally stable states. This is because either state is locally stable, small perturbation 
with respect to the steady state dies out in time. 


In comparison to the above model, it is easy to show that the potential function proposed 
in ref. [29] is also bounded and dependent on the total vehicle density N. It has a unique 
global minimum, which corresponds to the FT state at low density and HCT state at high 
density. However, the location of the minimum changes continuously with iV, so it does not 
present any sudden jump, neither the coexistence of different phases. 

Generally speaking, when it is necessary, one can introduce more local minima into the 
potential function to model other local stable states. We shall not pursue such complication 
any further, but turn to discuss the possibility to obtain non-stationary solution in this 
one-dimensional model. At a hrst thought, one might conjecture that an oscillating state 
corresponds to some periodic transitions between different states while possibly passing 
through some unstable local maxima. However, it is not difficult to show that the above 
conjecture is not possible: in a one-dimensional case, the system starts to evolve from any 
state, the nature of the EoM forces it always to evolve in a given direction in the conhguration 
space without ever turning back, and the evolution will only stop when it encounters a stable 
state. So unless there is a periodic external perturbation to constantly push the system away 
from a given stable state, the motion can not be periodic. In our model, since the stochastic 
noise does not possess any periodicity by dehnition, even if it were taken into account, it 
might not have served as the cause of possible periodic solution. Therefore one arrives to 
the conclusion that in the one-dimensional case, there is no periodic oscillatory solution. 

In what follows, we calculate the fundamental diagram obtained from this one-dimensional 
model, and present it in Figj^a). In our calculation, for simplicity, we assume Vi = 0, ^2 = 1. 
The details of the formulae can be found in the Appendix. The calculated fundamental 
diagram shows clearly an inverse-A shape. This is owing to the coexistence of the metastable 
states (dashed or dotted curves) and the stable states (solid curve), as discussed in the text 
concerning Figj^ The maximum flow of the FT state (the tip of the inverse-A) corresponds 
to a metastable FT state, while the corresponding stable steady solution is represented by a 
stable HCT state lying vertically below it: perturbations big enough may trigger the system 
to jump from the tip of the inverse-A to the HCT state. The jump is due to the difference 
of the magnitudes of the corresponding traffic flow of the two states, which is absent in the 
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FIG. 3. (Color Online) Fundamental diagram of the present traffic model: (a) Fundamental 
diagram of the one-dimensional case where the inverse-A feature is reproduced owing to the jump 
between the metastable and stable states. Here one uses solid lines to indicate stable states and 
dashed or dotted lines to indicate metastable states; (b) Fundamental diagram of the first scenario 
in the two dimensional case. The scattering of the flow-concentration relationship is due to the 
stable limit circle solution. When one removes the curl field B = V x and the auxiliary 
potential , the resulting fundamental diagram reduces to that of the the inverse-A shown in 
(a); (c) Fundamental diagram of the second scenario in the two dimensional case. The scattering 
of the flow-concentration relationship is due to a finite region of metastable states. Again, if 
one removes those metastable states, the resulting fundamental diagram reduces to that of the 
one-dimensional case shown in (a). 


model proposed in [29]. However, since one is only considering the deterministic EoM, the 
current approach does not contain a mechanism within itself to provide such perturbation. 
We note that such sudden jump of physical quantities appears naturally in our approach, 
in the sense that it is not originated from any discontinuity in the parametrization of the 
model, but from its nonlinearity. However, another important feature of the fundamental 
diagram, namely, the observed scattering of the data of the flow-concentration relation in 
the congestion region (as shown in Figj^b) and Figj^c)), is not reproduced here. Intu¬ 
itively, such feature may partly be related to the stochastic nature of the traffic flow m- 
Interestingly, as discussed below, it may also be attributed to either a stable limit circle 
solution or to a finite range of metastable solution in the parameter space. However, in the 
one-dimensional case, neither of the two solutions is attainable, they only exist in the higher 
dimensional case, as to be discussed in the following section. 

To close this section, we give a few more comments on the connection between the “cusp 
catastrophe” model and the Landau’s theory of phase transition. According to Landau [13] , 
the stability of a thermodynamic system at fixed temperature and volume is governed by its 
free energy, which can be derived from the partition function following standard procedure 
of statistical physics. In this phenomenological approach, the free energy of the system 
can be expressed as a function of order parameter T. A metastable state corresponds to 
a local minimum of the free energy, while a macroscopically stable state is encountered 
at the global minimum of the free energy. If the free energy has only one minimum, the 
thermodynamical system possesses only one phase. At the moment when a second phase 
begins to emerge, the extrema of the free energy become threefold, namely, two minima and 
one maximum. A phase transition takes place when the free energies of the two phases are 
equal. The phase transition is said to be of second order if the two phases coincide at a 
critical point, otherwise the phase transition is of first order. In the catastrophe model, it is 
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FIG. 4. (Color Online) Extrema of the corresponding “cusp catastrophe” model in its parameter 
space. Solid line corresponds to global minima, while dashed and dotted lines indicate local maxima 
and minima respectively. For instance, at = 0.35, the extrema is threefold, its potential function 
is shown in Fig{^b); while at = 0.6, the extremum is onefold, its potential function is shown in 
Fig^d). 

assumed that the physical state is attained at the global minimum of the potential function. 
Apart from the fact the potential function is usually an a priori ingredient of the model, the 
mathematics of the two approaches are strikingly similar. In the case of “cusp catastrophe” 
model, the potential function is threefold in some specihc region of the parameter space and 
onefold otherwise. To illustrate the above discussions clearly, the extrema of the above one 
dimensional model is shown in Figj^ where the positions of the extrema are determined by 
steady states of the equation of motion Eq.(|^. At A^ = 0.35 for instance, two of the three 
extrema correspond to stable solutions while the other one is linear unstable. In our model, 
the potential function is derived in terms of transition coefficients by Eq.([^, therefore it 
provides a mesoscopic interpretation for the catastrophe model stemmed from the dynamics 
of the traffic system. 


IV. HIGHER DIMENSIONAL CASE AND TEMPORAL OSCILLATORY 

SOLUTION 


The above discussions on scalar potential function can be generalized to the higher¬ 
dimensional case. The motivation to study the problem in higher dimensional parameter 
space is to investigate the non-stationary solution of traffic system. To increase the dimension 
of the parameter space, one notes that the set of equations Eqs.(|^ introduced above can 
equivalently be seen as a system with two degrees of freedom if no constraint is applied. 
Following the same arguments of the previous section, one can write down the EoM in 
terms of the gradient of a potential function U and dehne it as a vector held E as follows 


drii 

dt 

dn2 

dt 


El 

E2 


E = {Ei,E2) 


( 10 ) 


with the vector held E dehned as the gradient of the potential: 


E = -Vt/ 


(11) 
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FIG. 5. (Color Online) Schematic potential function of the two-dimensional case and possible so¬ 
lutions: Left: An illustrative two-dimensional potential featured by more than one local minimum: 
the minimum at the center of the potential represents the FT state, other minima correspond 
congestion states; Middle: Schematic temporal evolution in the configuration space. The red and 
black curves indicate the temporal evolution of the system against small deviations from FT and 
congestion states respectively, in both cases, the system quickly converges to the nearest stable 
state. The green circle indicates a possible scenario of a limit circle solution in the two-dimensional 
parameter space; Right: Semi-stable limit circle solutions of Eqs.(12). 


It is noted that a given state of the system corresponds to one point in the conhguration 
space; the system evolution is completely determined by the vector held alone, which draws 
a curve starting from the point in the conhguration space corresponding to the initial state. 
The evolution of the system can be viewed equivalent as that of a mechanical system ac¬ 
cording to Hamilton’s equation in classical mechanics. In Figj^ we depict in the left panel, 
an illustrative two-dimensional potential and in the middle, how the system can evolve in 
time in the conhguration space. The schematic potential has two minima. The hrst one 
is at the center of the potential, let us assume that it corresponds to the FT state at low 
vehicle density; the second minimum locates in the basin that circles around the center, 
which represents the congestion state; at the outskirt of the potential, the gradients are all 
pointing inwards, so that the system will not diverge. Apart from its complication, such a 
two-dimensional potential does not introduce new physics. As shown in the middle panel of 
Figj^ for small deviations from either FT or congestion state, the system quickly converges 
to the corresponding stable state. 

However, the green circular trajectory in the middle panel of Figj^ indicates a novel 
possibility. Unlike the one-dimensional version of the model, now it is possible for the 
system evolution to trace out a closed curve, therefore to form the so called limit circle in 
the parameter space as shown by the circular trajectory in the plot. Once the system evolves 
back to its initial condition, the EoM implies that the system will repeat the trajectory and 
carry out a periodic motion in the parameter space. Theorem on the existence of the limit 
circle is described by the Pioncare-Bendixson theorem [12] : roughly speaking, when no stable 
state exits, the system evolution must be featured by periodic orbit. Unfortunately, traffic 
system possesses at least the FT as a stable state, and therefore a general theory on the limit 
circle of such a system is yet unknown. In fact, for the specific form of the potential function 
shown in the left panel of Figj^ it is easy to show that the above periodic evolution will not 
occur at all when the vector field is solely determined by the gradient of the potential in 
Eq.(|g). This is because the potential in Eq.(pT| is conservative by definition, and therefore 
the corresponding field is irrotational. Nonetheless, an oscillatory solution can be explicitly 
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obtained if one introduces some rotation into the EoM as in the following equations: 



Tj j 2 I 2 Q 

= rii = B 2 = y nf + n 2 cos & 

9 = atan{n2/ni) 

Since the above vector held is solenoidal (divergenceless), it can be written as follows 


( 12 ) 


B = {B,,B2) 

B = V X A (13) 

where A is a vector potential in corryarison to the scalar potential U introduced before. 
Due to the fact that the vector held B is solenoidal, according to Pioncare-Bendixson theo¬ 
rem, it possesses limit circle solution which go round and round in the conhguration space 
perpetually. For such a simple case, the corresponding analytic solution is shown in right 
panel of Figj^ and is known as semi-stable limit circle. In a realistic case, an additional 
scalar potential may also exist. The behavior of the solution depend on the relative strength 
between the scalar and vector potential. If one assumes that the scalar potential is strong 
enough and its stable minimum locates at the center of the circular evolution, any small 
perturbation may cause the trajectory to eventually sink into the stable minimum. Such 
solution is known as unstable limit circle. Inhuenced by the shape of and relative strength 
between the vector held and scalar held, the circular solution may also be stable, so that all 
nearby trajectories spiral towards it. In fact, we understand that most empirically observed 
periodic oscillations are stable against small perturbations, and therefore must correspond 
to stable limit circle as discussed in ref. [20]. It is observed that one also encounters cases con¬ 
cerning the transition between periodic oscillatory state and diherent steady stable states. 
For instance, by adjusting the vehicle how at the entrance of the ramp, one may force the 
traffic phase to transfer from a stable FT into either an OST or a HCT [S] [201 E2|- Such 
phenomena can be interpreted as the transition between stable limit circle and stable critical 
points. Additionally, it is known in vector calculus [33| that a general vector held can always 
be decomposed into a sum of an irrotational held and a solenoidal held, which ensures that 
the above discussions is generally valid since the FoM can be always written in the following 
form 


dni 

dt 

dn2 

dt 


Ei + B, 
E2 B2 


(14) 


where the vector helds E and B are determined by the scalar (Fq.(ll)) and vector (Fq.(13)) 
potentials respectively. Roughly speaking, the scalar potential dehnes the locations of possi¬ 
ble sinks, while the vector potential creates possible limit circles, the interplay between the 
two results in diherent types of limit circle solutions accompanied by steady solution. 
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In what follows, we construct two different scenarios, where the model consists of an 
irrotational held E and a solenoidal held 5 in a two-dimensional case. In the hrst scenario, 
the model is aiming at reproducing a traffic theory consisting of FT, HCT and OST as studied 
in |22]. It is done by properly choosing the form of the scalar and the vector potential to 
form a stable limit circle for the OST state. In the second scenario, we modihed our model 
in order to incorporate the “synchronized how” state and the corresponding F—)-S transition 
introduced in the three phase traffic theory. There, the vector potential is essential for the 
oscillatory behavior before the synchronized how state is attained. 


The first scenario: an implementation for FT, HCT and OST states 

In this scenario, the potential function consists of a scalar potential and a vector potential. 
We relegate the detailed form of the potential functions to the Appendix, here we only 
outlines the rule of thumb for the choice of the potentials: (1) The scalar potential has three 
minima, which give rise to FT, HCT and OST respectively; (2) At low/high vehicle density, 
the minimum of the scalar potential corresponding to the FT/HCT dominates, so that the 
resulting potential is essentially the same as in the one dimensional case discussed in section 
HI; (3) At intermediate vehicle density, an oscillatory solution is found which is caused by 
the interplay between the scalar and the vector potential where the latter is the strongest 
in the vicinity of the corresponding minimum of the scalar potential. The latter is tuned 
to be strong enough to attract the oscillatory solution but not too strong to destroy the 
stability of the limit circle. (4) The range of the above stable limit circle shall reproduce 
the scattered feature of the observed data. 



FIG. 6. (Color Online) The two-dimensional scalar potential functions for a three-phase traffic 
theory: Left: The scalar potential at low vehicle density N = 0.2, it contains the global minimum 
corresponding to the stable FT state; Middle: The scalar potential at high vehicle density N = 0.8, 
it contains the global minimum corresponding to the stable HCT state; Right: The scalar potential 
at intermediate vehicle density N = 0.6 featuring two local minima, their interplay with the vector 
field B forms a stable limit circle, the OST state. The analytic expression of the potential function 
is given in the Appendix. 

In Figj^ we show the scalar potentials of the model for three different traffic states cor¬ 
responding different vehicle densities. To show the above model indeed describes the above 
different traffic states and their transitions, we show the corresponding temporal evolutions 
of the system in Figj^ The results are obtained by numerical calculations. To simplify 
the calculations, we assume that the velocities corresponding to the two free parameters 
rii and n 2 are almost degenerate, namely, ui ~ U 2 — 0 , while the third velocity U 3 = 1 . 
The hrst column of Figj^ depicts one of the parameters, the ratio between the occupation 
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density ni and the total density iV, as a fnnction of time. One notes that other physical 
qnantities (snch as overall vehicle density and flow) snbseqnently follow the same pattern. 
In the second colnmn we show how both parameters, ni/N and n 2 /N, evolve in the con- 
hgnration space. In the plots, different initial conditions are considered. It is seen that the 
characteristics of all three phases are well reproduced. In the cases of FT and HOT, the 
system always shows non-oscillatory relaxation to the corresponding stable state. On the 
other hand, at intermediate vehicle density, different phases coexist: some initial conditions 
can involve into OST where the system oscillates and approaches the limit circle, while the 
others with bigger initial flow are instead drawn into the FT state. The oscillatory solution 
in the third column of Figj^ is identihed as OST since the flow oscillates between a value 
(close to that) of a HOT state and an intermediate value. It is emphasized since the present 
model is introduced in a very generic context, the resulting limit circle solution is held in a 
general context which does not depend on specihc choice of parameters. For instance, TSG 
can also be easily incorporated into the model if one further introduces another additional 
rotation involving the minimum of FT state. 
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FIG. 7. (Color Online) Temporal evolution of the system for the first scenario of the two- 
dimensional case. From the left to the right are the three scenarios shown in FigJ^for low vehicle 
density, high vehicle density and intermediate vehicle density respectively. First column: Tempo¬ 
ral evolution of one of the parameters of the system. For the FT and the HOT case, the system 
converges to the corresponding stable state, for the OST case, the system converges to a periodic 
motion; Second column: Evolution of the system in the configuration space. For the limit circle 
solution, we note that the system never converges to any fixed state. 


The resulting fundamental diagram is shown in Figj^b). It is pointed out that due to 
our intended choice of parameters, the curve of FT and HCT states lie exactly at the same 
place as the one-dimensional version of the model, so that the difference between the two 
models becomes more transparent. In the present case, more phases are introduced into the 
fundamental diagram owing to the complication of the potential functions. The appearance 
of the grey band in the fundamental diagram is due to the oscillatory solution of OST, which 
occupies a hnite area in the diagram instead of a thin curve. Therefore it may account for the 
scattered feature of the observed data. It is also noted that stability of FT and HCT phase 
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is modified in some regions, although the positions of the curves remain unchanged. This 
is because the existence of a the HCT phase affects the stability of the others. As a result, 
part of the stable HCT phase (marked in solid red curve) in Figj^a) turns into metastable 
phase (marked in dashed red curve) in Figj^b); and the part of HCT overlapping with the 
OST phase simply disappears due to the strong influence of the vector potential. 


The second scenario: an implementation of synchronized flow state 

For this scenario, the synchronized flow state possesses the following characteristics p.lM]: 
(1) According to the three phase traffic theory, the synchronized flow state is related to 
metastable states. (2) The collection of metastable states accounts for the scattering data 
points in a region on the fundamental diagram where the traffic breaks down. As a result, 
there are infinite number of metastable states in the parameter space. In our model this is 
implemented by introducing a circle of metastable states in the parameter space, which are 
characterized by the same depth of the scalar potential function. (3) In order that there is a 
high probability for the transition F—)-S, in the region where free flow is a metastable state, 
the congested state is either an unstable state, or it is a metastable state but geometrically 
difficult to be reached directly from the free flow state in the parameter space. In both 
cases, the chance to observe a direct F—)-J transition becomes impossible or insignificant. In 
our model, the first possibility is adopted. (4) The vector potential introduces oscillatory 
behavior in the F—)-S transition process. 

Again, the specific forms of the potential functions are relegated to the Appendix, we only 
note that the form of the potential is chosen so that the inverse-A shape of the fundamental 
diagram (Figj^c)) remains at the same position as in previous cases. For simplicity, the 
region of metastable states is taken to be a closed circle in the parameter space of the 
model, as illustrated in Figj^a). In Figj^b) and (c), we show the temporal evolution of 
the system for different initial conditions in the synchronized flow phase. As shown in 
Figj^b), depending on the initial condition, it is shown that the resulting state can pick 
any one from the family of metastable states shown in Figj^a). If some small fluctuation 
is introduced, the system will be perturbed but evolves and eventually falls into another 
metastable state. In Figj^c), we show the corresponding fundamental diagram. In the plot, 
the region featuring the scattering of the data is due to the region of metastable states. In 
terms of the fundamental diagram, it seems that both scenarios capture the feature of the 
scattered data. Figgc ) also shows clearly that it is impossible for the system to evolve 
directly from a free flow state to the congested state. This is because for a given vehicle 
concentration, namely, a vertical line, will never intersect with both the curve representing 
the free flow state and that of the congested state. 


V. CONCLUDING REMARKS 

In this work, a mesoscopic transport model is employed to study different traffic phases 
and their transitions. Our approach is based on a mesoscopic gas-kinetic viewpoint of the 
traffic system in the same sense as proposed by Prigogine and Andrews [35]. By defining 
the potential function in the catastrophe model, the present work establishes an intuitive 
explanation for different states in traffic flow and their connection to the stability of the 
corresponding EoM. Our approach bears close connection to the catastrophe model and 
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FIG. 8. (Color Online) Metastable region of the potential function and temporal evolution of the 
system for the second scenario of the two dimensional case. Left: the minima of the potential 
function in the parameter space; the synchronized flow correspond to the circular region contains 
an infinite number of states, while other point-like stable states such as FT and HCT are also 
marked on the plot for comparison; Middle: Temporal evolution of one of the parameters of the 
system. Depending on the initial conditions, the system oscillates and converges to any point in 
the metastable region where III and III indicate three different initial conditions; Right: The same 
as the middle plot, the evolution of the system is shown in the configuration space. 


Landau’s theory of phase transition. In the one-dimensional case, by adopting an appropriate 
form of the potential function, the model is able to reproduce stable steady traffic phases, 
including FT and HCT. It reproduces the inverse-A feature in the fundamental diagram. 
When generalized to the higher-dimensional case, it is shown that the model provides a 
mechanism to describe the temporal oscillation in conhguration space. In our interpretation, 
the OST state is attributed to the formation of a limit circle; the synchronized flow state 
is attributed to a finite region of metastable states of the potential function. Therefore, 
elaborating such a dynamical theory not only provides a possible microscopic origin behind 
the potential function, but also makes it feasible to intuitively deal with phenomena such as 
stationary as well as non-st at ionary stable states. The introduction of the potential function 
helps not only to quickly identify the locations of metastable states as well as oscillatory 
states, but also to quantitatively investigate the stability of such states. 

As mentioned above, the topic of traffic flow states and their transitions have been ex¬ 
plored by many authors. Most of the works employ either macroscopic hydrodynamics or 
microscopic car-following models. Such approaches are usually based on detailed observa¬ 
tions including traffic data and/or realistic highway and vehicle parameters. Due to such 
sophistication, heavy numerical simulations are usually required where a thorough scanning 
of the model’s parameter space as well as the initial conditions become increasingly diffi¬ 
cult. In addition, specific details in the construction of the model lead to delicate differences 
in the interpretation of numerical results, especially when the resulting solutions are not 
unambiguously identified in the observed data. On the other hand, in the present study, 
we aim to develop an approach from a different angle. Instead of focusing on specific de¬ 
tails in modeling the traffic system, we concentrate on how the resulting traffic states are 
connected to the concept of instability and phase transition. Inspired by the catastrophe 
model, all possible stable steady traffic states are associated with the local minima of the 
scalar potential, while vector potential gives rise to possible oscillatory traffic states. The 
simplihed mesoscopic gas-kinetic theory introduced in this work is to show the underlying 
physics more transparently. In particular, most of our discussions are from a phenomeno- 
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logical viewpoint. In this context, many findings of the current study are not dependent 
on any specific traffic scenario, neither on the particular external source in the EoM nor 
specific boundary conditions. It is an impending task to further apply the present model to 
the description of the data. 

In the above higher-dimensional models, the wide scattering of traffic data in the fun¬ 
damental diagram is either interpreted in terms of the OST phase corresponding to a sta¬ 
ble limit circle, or in terms of the synchronized flow phase corresponding to a region of 
metastable states. For the first scenario, the OST phase is usually obtained by introducing 
pulse-type perturbations |8], |T9l EHl l22] into an initially homogeneous metastable state. In 
our approach, coexisting solutions are obtained by choosing the initial conditions (as shown 
in the right column of Figj^ once the vector potential is properly given. As we understand, 
both approaches effectively kick the system into a certain region of its configuration space, 
therefore bear the same physical content. However, the cause of the limit circle is attributed 
to the curl of the potential function rather than to the inhomogeneity. In our model, once 
the shape of potential functions are determined, it is relatively easy to find out what specific 
initial condition may lead to limit circle solution. It is understood in [19] that oscillatory 
traffic state far away from the ramps are transient process of decaying oscillations. In terms 
of our interpretation, this happens when the depth of the minimum of the corresponding 
scalar potential is so big that the limit circle nearby becomes unstable. Therefore such phe¬ 
nomenon can be naturally described within the same framework. For the second scenario, 
the synchronized flow phase is obtained by introducing a region of metastable states, in the 
same sense as implemented in Ref.[2S|- According to the three phase traffic theory EHIEZ], 
the highway capacity is not a well-defined physical quantity. To have a region of metastable 
states implies that different values of highway capacity for different highway capacity, there¬ 
fore it is consistent with the above assertion. In addition, as shown in Figj^c), our approach 
adopts an adaquate set of parameters that a direct transiton from free flow to congested 
flow becomes impossible. Comparing to Figj^b), this can be achieved by decreasing the 
range of metastable congested flow. This can be done without much difficulty since the 
depth of scalar potential is easily manipulated in our model. In fact, by modifying the 
depth of corresponding metastable states, one can further adjust the probability of traffic 
break down as well as of the transition rates between different traffic states in accordance 
with experimental observations. 

As discussed before, there are two possible sources of instability. This work only involves 
one of them, namely, the stability of the deterministic FoM and its relation to the traffic 
flow. However, as traffic flow is stochastic in its nature, and our model itself involves a 
mathematical description of uncertainty, it is intuitive to attribute the observed randomness 
in the data to the stochastic noises. Such randomness, as we understand, can naturally 
provide a perturbation that takes place in reality to trigger the transition between different 
traffic phases. Besides, with the presence of noises, phenomenon such as hysteresis effect 
|2S| is expected to be reproduced. In fact, the problem of stability of SDE has always been 
a topic of increasing attention in its own right [SSHHj • It is meaningful to further construct 
a more realistic version of the model with its parameters determined by the observed char¬ 
acteristics of the traffic system. We note that a preliminary attempt was made recently BZ], 
The mesoscopic model for traffic flow has also been used to derive the corresponding macro¬ 
scopic hydrodynamical equations |13H15|. The same procedure can be readily borrowed to 
derive the corresponding hydrodynamical equations and carry out studies of the traffic flow 
system from a macroscopic viewpoint as the latter is closely associated with experimental 
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observables. It is interesting to further investigate and compare the stability of the traffic 
ffow states in both approaches. Do the spatial-temporal properties of the obtained hydro- 
dynamical model reflect realistic characteristics of the traffic system? Moreover, are they 
related to those of the corresponding Boltzmann equation? These are worthy topics for 
further exploration. 
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VII. APPENDIX: THE POTENTIAL FUNCTIONS EMPLOYED IN THE WORK 

For the first scenario of the two-dimensional case, the scalar potentials shown in Figj^ 
read 


r/( 2 ) 

vv ’ n) 

(2) frh r^\ 

* \n'nJ 


Tt( 2) frn r^\ 

^ \n'n) 



\N'n)^ ^ \N'NJ 


,( 2 )^ frii n2 \ 

= dl 'G 2 


d = 1,2,3) 
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While the vector potential in 


Eq.(13) is defined as 
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where the scalar and vector potentials are defined in terms of the ratio between occupation 
density of the two low velocity states ni 77-2 and the total vehicle density A, points in the 
direction perpendicular to the parameter space. The parameters, the Gaussian function G 2 
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and Kernel fnnction K 2 are defined as follows 


with 


and 


df) = -0.01K:2(Ar,iV,0,0,0.4), 
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df) = -0.05K:2(A^, iV, 1,1, 0.6), 

= 700df\ 


Xi = Ui = 0, cTi = 0.05, 

X2 = y2 = c^s = 0.03, 

1 + N 

X3=y3 = -, 0-2 = 0.04, 

3 1 + iV 

Xr = yr = -1-, CTr = 0.5, 


20 10 


a = 0.1 


G 2 {x,y,Xo,yo,a) = exp 
V 27r 


= 5^(2-?)= 

0 


(2;-a;o)2 + (2/-2/o)^\ 

(18) 

2a2 J 

(19) 

'(x - 2 : 0 ) 2 + (?/- 2 / 0 ) 2 , hj 

(20) 

- 0 < g < 1 
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As discnssed in the text, the scalar potential which featnres three minima, corresponding 
FT, HCT and OST states. In onr calculations, for simplicity, the maximal value of N is 
scaled to be 1.0. The form of the minima is produced by a Gaussian function. On the one 
hand, the minimum Xi is related to the FT at low vehicle density, its position is hxed in 
the parameter space so that flow-concentration curve has a constant inclination for the FT 
phase as shown in Figj^b). As discussed in section III, when the overall vehicle density 
increases, the minimum of FT starts to become metastable and eventually disappears. This 
is controlled by the Kernel function K 2 defined above. Similarly, the position of minimum 
X 3 corresponding to HCT moves towards 0.5 (as {jii + n 2 )/N moves towards 1.0) as N 
increases, so that more and more vehicles transit to the low velocity state ni and n 2 which 
consequently causes the overall traffic flow to decrease. Since the two minima satisfy 2:3 > Xi 
in the region where FT and HCT coexist, the congested flow is always smaller than the free 
flow, and the difference corresponds to the jump as shown in Figj^ On the other hand, the 
additional local minimum X 2 is introduced here to form the OST phase, together with the 
localized rotational held = V x dehned in terms of the vector potential Since 
the Kernel function K 2 is centered at = 0.7, the strength of the vector held B is also a 
function of the total vehicle density iV, which is suppressed at lower and higher total vehicle 
densities. It is straightforward to verify that the vector held B forms closed circles around 
{xr,yr) with an ehective range of It is noted that even though the equation of state of 
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our two dimensional model is defined to be symmetric between ni and n 2 , it solution, as 
shown in Figj^ depending on the choice of initial conditions, does not necessarily possess 
such symmetry. To numerically evaluate the resulting traffic flow, one assumes for simplicity 
that the two low velocity states are almost degenerate, namely, ui ~ U 2 — 0 , while the higher 
velocity U 3 = 1 . 

The one-dimensional case can be obtained by removing the vector potential as well as 
the minimum corresponding to the OST state, and also considering only one unique low 
velocity state (perfect degeneracy) rii. Due to the following intentional choices of the pa¬ 
rameters, the one-dimensional model can be simply seen as a special case of degeneracy of 
its more sophisticated two-dimensional correspondence, and therefore corresponding flow- 
concentration relation for FT and HCT are exactly the same in the two plots in Figj^ To 
be specific, in the one dimensional case, the scalar potential is defined as 



( 22 ) 


where 


df = 0, 

42 ^ = 0, (23) 


and all other parameters remain unchanged. Therefore the two minima Xi and are 
interpreted as FT and HCT states, in exactly the same way as in the above two-dimensional 
case. 

For the second scenario of the two-dimensional case, the scalar potentials are essentially 
the same as those of the first scenario, an additional region of metastable states is introduced 
in the following: 
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While the vector potential in Eq.(13) is defined as 


1(2) = H(2)e, 
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where the parameters are as follows 

G2 - Xr^ + {y- VrY, x„, x„, aj 


F =1 

^ a -*- 


rfp) = -0.01K2{N, N, 0, 0, 0.4), 
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According to our mesoscopic flow, the traffic flow can be therefore obtained in terms of 
the expected value 

g = ^ TiiVi (28) 

i 

Thus the flow-concentration relation can be evaluated once the temporal evolution of Ui is 
obtained numerically. It is noted that the soluction of the one dimensional model can simply 
be obtained by substituting (ni -|- n 2 ) —)■ ni and n 3 —)■ n 2 in the two dimensional case. 
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